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Thanks to Wolstenholme m, the following three congruences have 
been know since 1862, for all primes p > 5: 

(!) _ J 1 ) = 1 ( mod P 3 )> 

(2) 1 + ^ + ^ + • • • H-—r = 0 (mod p 2 ), 

2 6 p — 1 

(3) 1 + h + h + '" + ( P -iy =° (modp) - 

Here, of course, | means the (multiplicative) inverse of k in the relevant 
sense: in 7L V , Z p 2 , etc, according to the context. More than 125 years 
later, Gardiner [1] showed the relation between these equivalences when 
the degree is pushed one level higher: 

Theorem 1. If p > 7 is prime, the following conditions are equivalent: 

(a) p is a Wolstenholme prime, meaning : f 2 ^ 1 ) = 1 (modp 4 ), 

(b) 1 + | + | + -- - + ^zt = 0 (mod p 3 ), 

(c) 1 + £ + p + • • • + ^zTj 2 = 0 (mod p 2 ), 

(d) p | B p _ 3 where Bdenotes the k th Bernoulli number. 

The term Wolstenholme prime was introduced by McIntosh in [10] . 
The only known Wolstenholme primes are 16843 and 2124679; see [$| 
for another equivalent condition. Nevertheless, Gardiner’s result has 
been extended one degree further. 

Theorem 2. Ifp > 7 is prime, the following conditions are equivalent: 

(a) ( 2p “ 1 ) = 1 (mod p 5 ), 

(b) l + 1 + i + ... + ^^O (mod p 4 ), 

(c) 1 + ^ + ^ + ... + ^^2 = 0 (mod p 3 ), 

(d) p 2 | B p 3_p2_ 2 . 

The above result is implicitly contained in Helou and Terjanian’s 
2008 paper [4], but somewhat scattered amongst a raft of other, often 
more substantial results. We will say more on this at the end of this 
note. Our main goal here is to highlight the result itself, and to provide 
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a unified, elementary and direct proof. One basic classical result we 
use freely throughout this note was proved by Leudesdorf [5]; see also 
[31 Chap. VIII.8.7] and 0: 


Lemma. If p> 7 is prime and k E N such that 2k < p — 1, then 

(4) ^ ii = ° (mod p 2 ), and 

l<i<p-l L 

(5) ^k=° (™dp). 

1 <i<p— 1 


Proof of Theorem [H (a) <=> (c). We first develop the binomial coeffi¬ 
cient ( 2 pli ) “downwards”: 

'2p - 1\ = (2p-_ 1)(2 p - 2)... (2 p-(p- 1)) 

p-l) 

2 p 


1-2. ..(p-l) 

Expanding the last line in Z p s gives us: 

1 

ijk 


p — 1 


(6) l-2p^^ + 4p 2 ^^-8p 3 £ 7^ + 16/ E ^jr v 


i i<j i<j<k i<j<k<l 

where here and below, unless otherwise stated, the summations are 
over variables in the range 1, ...,p — 1. Next we work “upwards”: 


2p- 1 
p — 1 


to obtain in Z p 5 : 


(1 + p)(2 +p) . . . ({p- 1) +p) 


1-2...(p-l) 


V 


= 1 + T 1 + X 


P 


1 + 


P 


P — 1 


(7) ^+pJ2f + p 2 J2h +p3 E fjk+p* E 




l 

ijkl 


m 


i<j u i<j<k i<j<k<l 

Multiply equation (J7D by 2 and add the product to equation 
order to eliminate the p term. Then divide both members by 3 to get: 

' 2 p s i -+:v £ £ -- ip £ ff k +'<#* £ pr 

7 Kj i<j<k J i<j<k<l J 


Concerning the last summand, notice that multiplying all the indices 
i.j. k, l by 2 leaves the sum ijg fi xed i 11 Therefore, since 

2 4 ^ 0 (mod p), this sum is equivalent to 0 (mod p). The second 







WOLSTENHOLME AGAIN 


3 


summand may be transformed by using 2 ^2 -4 = |) 2 — Y2 -j. After 

substitution and application of (J2]) to the square term we get: 

(T-'i 1 ) - 1 - p 2 £ ^ - ■ 2 p 3 £ ^ (>“ od p 5 )- 

x ^ 7 i i<j<k J 

Finally, concerning the last summand, notice that we have: 

(£j)+*£* 

i<j<k lJK \ i V V < 1 / V j 3 J i 1 

which is equivalent to 0 (mod p 2 ) by using the equivalences (J2j), (J3j) 
and (pi]). Therefore we have proved 

(s) ( 2 p p _ ^ = 1 - p 2 Z p ( mod p 5 )’ 

which figures in [10, p. 385]. 

(b) <(=)> (c). By using elementary identities we obtain: 

^Z ((p — i)i i 2 i 2 ) 

= -f>E^+p 2 E(^2 

z i x ' 

= -pJ2^ + p E ((p~i)i 2 + p) Zp 

z z x v 7 7 z 

= -?E^-p 2 E^+j> 3 E(p^. 

z z z x 7 

from which we easily conclude by using fl3D on the middle summand 
and ([5]) on the last summand as JA = JA (mod p). 

Equation (J8]) <(=)> (d). This last equivalence requires basic knowledge 
of Bernoulli numbers we recall from [5]. If 

p -1 

(9) S m (p):=Y,i m 

i =1 

then from [5] pg. 230, Theorem 1], 

m+ 1 / \ 

(10) S m (p) = ^ 


2et=e 


p 


1 1 
- + - 
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Importantly for us, the B t are 0 for odd integers i > 1. Our general 
method is to transform the summand in (JS]) into an equation of the 
form ([9]) by applying Euler’s theorem, 

i~ 2 = ^ (p3) ~ 2 (mod p 3 ), 


where (f) is Euler’s totient function. Working in Z p 3 and letting m : = 
p 3 — p 2 — 2 we get, since odd indexed Bernoulli numbers vanish, 




Smijp) = 


m +1 1 

E 1 

fct ! 


m 

i — 1 


P 1 B, 




= p Bp3 — p 2 —2 (mod p 3 ) 
which replaced in (JHJ) gives what is wanted: 

(^p_ i) = 1 _ P 3 B p 3- p 2_2 (mod p 5 ). 


□ 


McIntosh commented that there is probably only a finite number of 
primes verifying Theorem [2](a) and conjectured that there are none [ TUI 
bottom p. 387]. One natural question is: Can Theorem [2] be extended 
to the next degree? According to (U Lemma 3 and Cor. 5(1)] it seems 
the answer is no, since they obtain the following results: 


2p — l\ 16 

_ 1 ) = 1 - p 3 Bp3 — p 2 —2 + ^P 5 B P -3 - -P 5 5 p _ 5 (mod p 6 ), 


p - 1 


5^ \ = - P -B pi _p 2_2 + ^S p _3 - yfi P -5 (mod p 5 ), 
2=1 

1 p3 4 

X/ ?2 = P B p3-p2- 2 - y^p-3 + g A 3 -Bp-5 (mod p 4 ), 


2=1 


and so the last term in i ? p _5 does not coincide in any pair of expres¬ 
sions. Notice that reducing these three equivalences modulo p 5 ,p 4 ,p 3 
respectively establishes Theorem [2[ It is in this sense that Theorem [2] 
is contained in [4], A formula for (^y 1 ) modulo p 7 is given in |9j. For 
related results see [7j. 
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